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The outer crust structure and composition of a cold, non-accreting magnetar are studied. We 
model the outer crust to be made of fully equilibrated matter where ionized nuclei form a Coulomb 
crystal embedded in an electron gas. The main effects of the strong magnetic field are those of 
quantizing the electron motion in Landau levels and of modifying the nuclear single particle levels 
producing, on average, an increased binding of nucleons in nuclei present in the Coulomb lattice. 
The effect of an homogeneous and constant magnetic field on nuclear masses has been predicted by 
using a covariant density functional, in which induced currents and axial deformation due to the 
presence of a magnetic field that breaks time-reversal symmetry have been included self-consistently 
in the nucleon and meson equations of motion. Although not yet observed, for B > 10^®G both 
effects contribute to produce different compositions -odd-mass nuclei are frequently predicted- and 
to increase the neutron-drip pressure as compared to a typical neutron star. Specifically, in such a 
regime, the magnetic field effects on nuclei favor the appearance of heavier nuclei at low pressures. 
As B increases, such heavier nuclei are also preferred up to larger pressures. For the most extreme 
magnetic field considered, B = 10^®G, and for the studied models, the whole outer crust is almost 
made of loZrsz. 

PACS numbers: 21.10.Dr, 21.30.-x, 21.60.Jz, 26.60.Gj, 26.60.Kp 


I. INTRODUCTION 

Pulsars have typical surface magnetic fields of 10^^ G, 
some displaying up to 10^^ G. Larger magnetic fields on 
the surface of neutron stars have been observed or derived 
from observational data, with a maximum inferred value 
of 2.4 X 10^^ G in soft gamma-ray repeaters (SGRs) and 
anomalous X-ray pulsars (AXPs) [TH5] . Magnetic energy 
might be released in star-quakes producing short bursts 
of gamma-rays observed in SGRs and AXPs [6]. The 
variation in the luminosity detected in such events sug¬ 
gests B ~ G m- Stronger interior magnetic helds 
are thought to be present as suggested by various obser¬ 
vations [SHlQ]. 

Based on theoretical calculations, the possible exis¬ 
tence of magnetars displaying larger magnetic fields have 
been hypothesized by considering dynamo-like effects on 
the interior of very young neutron stars. Such effects 
might be enhanced by considering magnetic instabilities 
during the supernovas just before the neutron star is 
formed [nHi3]. According to the virial theorem and 
magneto-hydrodynamic simulations, the upper limit for 
the neutron-star interior magnetic helds is about 10^® G 

[mils]. 

Matter properties in the outer envelopes of a magne¬ 
tars are thought to be signihcantly modihed by strong 
magnetic helds m- Magnetic helds alter transport pro¬ 
cesses and therefore transport properties, such as thermal 
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and electrical conductivity [TMTO] . In addition, magnetic 
held stresses are thought to produce seismic modes that 
can be observed from Earth as quasi-periodic oscillations 
(QPOs) in the X-ray hux of giant hares from SGRs |20j . 

For magnetic helds B < 10^® G, the outermost layer 
of neutron star spans about seven orders of magnitude 
in density: from 10^gcm“®, the complete ionization 
density, to approximately 4 x lO^^gcm’ HU El], the 
neutron-drip density. In the simplest model, this outer 
crust is assumed to be made of nuclei and free electrons 
in their ground state at zero temperature |16] . Further¬ 
more, we consider the system embedded in a uniform 
magnetic held. Within this typical density range, the 
ionized nuclei hnd energetically favorable to arrange 
themselves in a Coulomb lattice [2T| . In the lowest range 
of density, the energy associated to the electron gas and 
to the Coulomb crystal does not play a relevant role 
in the determination of the nuclear species present in 
the crust: thus, the most stable isotopes of nickel and 
iron will most likely show up in such conditions. As 
the density increases, the lattice energy effects remain 
negligible, while the energy of the electron gas raises 
signihcantly as compared to the total energy of the 
system. Such energy is lowered decreasing the electron 
number in the gas through electronic capture processes. 
These reactions imply a progressive neutron-enrichment 
for the nuclei in the lattice. Meanwhile this mechanism 
cannot go on indehnitely because of the growing nuclear 
symmetry energy. Thus the outer crust is the result of 
a competition between the electronic energy, favoring 
neutron-rich nuclei, and the nuclear symmetry energy, 
favoring fairly symmetric ones. The outer crust ends 
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when the nuclei become unstable against neutron emis¬ 
sion because of their high neutron imbalance, becoming 
the inner crust where a free neutron gas is also present 
[m [23]. The ideal boundary between the two crusts is 
called neutron-drip transition, and the related dripping 
pressure Pdrip plays a fundamental part in order to 
calculate the outer crust spatial extension [M] Hi] . 


Nuclei present in the Coulomb lattice might be very 
exotic. Rare Ion Beam Facilities worldwide aim at ex¬ 
tending current mass measurements up to extreme values 
of neutron and proton asymmetries. However, accurate 
nuclear mass models predict the appearance of nuclei in 
the crust that have not been measured yet. The neu¬ 
tron star outer crust composition has been studied in 
various works, both in the absence [23H2S] and presence 
of magnetic fields [H 123I11S[3I]- Nevertheless, none of 
these previous analysis accounted for the effect of the 
external magnetic field on the nuclear binding energy of 
the nuclei present in the Coulomb lattice in a fully self- 
consistent manner. It will become evident in Sec. lIIIl how 
this dependence can change the ground state properties 
of matter in the crust. Specifically, we will rely on two 
commonly used covariant energy density functionals that 
have been successful in the global description of nuclear 
masses, charge radii, deformations and nuclear collective 
excitations. One of the models is based on the non-linear 
Walecka model [5^H51] and the other, more modern, on 
an effective Lagrangian with density dependent meson- 
nucleon couplings [35] . 

The aim of this study is two-fold. First, to determine 
the composition of a cold non-accreting and strongly 
magnetized outer crust {B > 10^"* G) focusing on the dif¬ 
ferences on the composition and structure of the outer 
crust obtained when the effects of the magnetic field on 
nuclei are included or neglected. And second, to ascer¬ 
tain if there are clear and model independent signatures 
of these differences, and equally importantly, at what 
magnetic field magnitudes they become relevant. We 
will present our predictions up to the theoretically-based 
maximum magnetic field of H ~ 10^® G, that guarantees 
mechanical stability to the star [111131131]. 

The article is organized as follows. In Sec. [TT] the for¬ 
malism is briefly discussed. In Sec. Ill we present the 
results for the structure and composition of the outer 
crust of a magnetar under the effect of strong magnetic 
fields. Our conclusions are laid in the last section. 


II. OUTER CRUST IN STRONG MAGNETIC 
FIELDS 


completely ionized atoms (with proton number Z, neu¬ 
tron number N and baryon number A) and a uniform 
Fermi gas of relativistic electrons [M m]. Through¬ 
out the outer crust a continuous value of the pressure 
is required in order to ensure hydrostatic equilibrium; 
the matter density can suffer from discontinuities due to 
the change of the nuclei present in the Coulomb lattice. 
Hence, the Gibbs free energy per baryon g{A, Z\ P) at 
a constant pressure and zero temperature should be the 
thermodynamic potential to be optimized rather than the 
energy per particle e{A, Z] p) at constant matter density. 
At zero temperature 


9{A,Z-,P) 


E[A,Z-P)+PV 

A 


e{A, Z]P) -\ -(I) 

n 


where V is the volume occupied by a unit cell of the 
Coulomb lattice, n = A/V is the baryon density in such a 
cell and e{A, Z; P) = E{A, Z; P)/A is the corresponding 
energy per nucleon. The minimum of g{A,Z]P) at a 
fixed pressure determines the couple A, Z that provides 
the most stable configuration. The energy per baryon, 
£{A,Z]P), is composed of three terms. 


e(A, Z; P) = e„(A, Z) + ee{A, Z; P) + ei{A, Z; P) , (2) 


the nuclear, electronic and lattice energy terms per 
baryon, respectively. The nuclear term accounts in our 
model for the nuclear mass M{A, Z), 

M{A,Z) Zrup + {A - Z)mn BE{A,Z) 
e„(A,Z)=—J-^ 

(3) 

where rrip and are respectively the proton and neu¬ 
tron rest mass, BE{A,Z) is the binding energy and we 
use, hereafter, natural units (fi. = c = 1). The nuclear en¬ 
ergy does not depend on the density in a unit-cell since 
it is much smaller than the nuclear saturation density 
throughout the outer crust. Within our description, in¬ 
dividual nuclei do not contribute to the total pressure 
which is composed of electronic and lattice contributions, 
respectively Pe and Pi. It is customary to relate, via 
charge neutrality, the electronic Fermi momentum ppe 
with the baryon average density n and define an average 
baryonic Fermi momentum pp ss follows |28j . 

PFe = ^ =(^jy PF (4) 

where rig is the electron density. Such an explicit rela¬ 
tion will be useful analyzing the electronic and lattice 
energies. 


We will highlight in this section the main features of 
our model for the neutron star outer crust, addressing 
the reader to 1131137] and references therein for further 
details. 

The matter in the outer crust of a cold (T = 0 K) non¬ 
accreting neutron star consists of a Goulomb lattice of 


A. Effects of the magnetic field on the electron gas 

Extreme magnetic fields affect the electron gas energy 
and are able to change the outer crust structure and com¬ 
position [13] . The existence of the magnetic field requires 
a different treatment and, thus, a new analytical form of 
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the rie and g{A,Z\P) with respect to the B = 0 case 
[HES]. In the plane orthogonal to the magnetic field 
direction -we consider a uniform magnetic field B = Bz, 
the electron motion is quantized into discrete Landau lev¬ 
els [3H]. Assuming a relativistic electron Fermi gas em¬ 
bedded in an external uniform magnetic field, the Landau 
energy levels E{i/,pz) can be written as follows [^ . 

E\„,p,)=pl + ml{l + 2uB,) = (5) 


where Pz is the electronic momentum along the mag¬ 
netic field direction, me is the electron rest mass, E is 
the electron energy, v is & non-negative quantum num¬ 
ber and Bi, is the external magnetic field B defined in 
units of the critical magnetic field Be- The critical field 
is defined as the magnetic field at which the electron cy¬ 
clotron energy equals the electron rest mass energy. That 
is. Be = m1/e « 4.41 • 10^^ G. 

Comparing Eq. (§ with the usual relativistic energy- 
momentum relation, a third additional energy term ap¬ 
pears due to the electron interaction with the magnetic 
field. This interaction energy is proportional to the quan¬ 
tum number v, and cannot exceed the electron chemical 
potential /ig. Hence the maximum number of Landau 
levels J-'max, related to the highest value of the interac¬ 
tion energy allowed between electrons and the external 
magnetic field, is evaluated setting E(venax,Pz = 0) = /ig 
in Eq. ([^. This leads to the expression, 


^max 


1 



( 6 ) 


Erom this equation, one sees that t'max is inversely pro¬ 
portional to the magnetic field. The latter is defined as 
strongly quantizing if only the lowest level is filled. In 
the general case, the maximum electron momentum avail¬ 
able, that is the Fermi momentum p-pe, can be computed 
setting E{iy,pz = ppe) = Me in Eq. ([^ for different iz 
values as 


[pFe{l^)f + Wg(I -f 2vB^) = pI 0<v< r-max (7) 


It is customary to define the adimensional Fermi mo¬ 
mentum Xeiv) and the adimensional Fermi energy 7 g as 
Xe(y) = ppe{v)lme and 7 e = Me/nie, respectively. The 
electronic density rig, energy eg, and pressure Pg can be 
then calculated as M- 


He 


Bi,ml 

47r^ 


g,^Xe{v) 

v=0 



V = Q 


( 8 ) 


7^ 4 ^tnax 

Eg = ^ 5 ,,(I -f 2vB^) r+ 

1^=0 


27r^ 




- 2vB^ 


(9) 


Pg n^Pe 


Bi,ml 

2'k'^ 


^ 5^(1 -f 2vB^,) r_ 


;y=0 


' Xe{v) 

- \/l -f 2vBi, - 


( 10 ) 


where, 


T±{x) = ^x^/ 1 -|- ± ^ In (x -I- \/l -I- a;2) (11) 

Unlike the B = 0 case, these functions cannot be stud¬ 
ied analytically, with the exception of the strongly quan¬ 
tizing magnetic field case. A detailed analysis can be 
found in [25]. Finally, we note that we neglect the small 
electron exchange corrections |29j . 


B. Effects of the magnetic field on the binding 
energy of the nucleus 

Theoretical extrapolations on nuclear masses are re¬ 
quired in order to describe the equation of state of the 
outer crust. Calculations of the needed nuclear binding 
energies employed in this work (both in absence and in 
presence of magnetic field) are based on the relativistic 
mean-field effective interactions NL3 [22| and DD-ME2 
[35] . As mentioned, the former model corresponds to a 
non-linear Walecka model while the latter is based on 
an effective Lagrangian with density dependent meson- 
nucleon couplings. Differences between the models will 
allow us to assess -to some extent- the model depen¬ 
dence in our results. Within the original and subsequent 
works, NL3 and DD-ME2 have been shown to be accu¬ 
rate in the description of experimental data on binding 
energies, charge radii, quadrupole deformations and the 
excitation energy of nuclear Giant Resonances [31]. 

The effects of the magnetic field on the binding en¬ 
ergy of nuclei have been taken into account following the 
work of Ref. m- For completeness, we highlight here the 
main features and address the reader to this work and 
references therein for further details. We have adopted 
a covariant density functional based on an effective La¬ 
grangian with nucleons and mesons as the effective de¬ 
grees of freedom [40] [41] : 

/L — H” ^int + Cbo + ^bm', (12) 

where Cn refers to the Lagrangian of the free nu¬ 
cleon, is the Lagrangian of the free meson fields 
and the electromagnetic field generated by the protons 
and Lint is the Lagrangian describing the interactions. 
These three terms compose the standard relativistic La¬ 
grangian. Throughout this work, the parameter sets 
NL3 [33] and DD-ME2 [3S] will be employed. These 
models differ on the form of Lint- Specifically, for the 
case of NL3, 


Lint = gatpertp - ^ 520 -^ “ ^ 53 ^^ 

- ePpjf,A^ip , 
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and for the case of DD-ME2, 

- - c'tp'y^A'^'ilj , 

where e is the electric charge for protons, '0 denotes the 
Dirac spinor, 7 ^ the Dirac matrices, and t the Pauli 
matrices in isospin space. The meson-nucleon vertexes 
are denoted by for z = cr, w and p; scalar-isoscalar, 
vector-isoscalar and vector-isovector helds, respectively. 
For the case of NL3, gi are constants and two non-linear 
terms have been introduced in the a field when com¬ 
pared to the case of DD-ME2. For the case of DD-ME2, 
the coupling constants are assumed to depend on the 
baryon density [35]. The adopted ansatz for the density 
dependence has been guided by more fundamental Dirac- 
Brueckner-Harthree-Fock calculations in infinite nuclear 
matter. The total number of adjusted parameters to 
some selected experimental data are six for NL3 and eight 
for DD-ME2. 

In addition, there are two terms corresponding to the 
interaction of the nuclear system with an external mag¬ 
netic field. The coupling of the proton orbital motion 
with the external magnetic field, Cbo = — 
and the coupling of protons and neutrons intrinsic dipole 
magnetic moments with the external magnetic field [42j 

Cbm = where , 

external field strength tensor, = 

§[ 7 Ati 7 i/]; MAf = ehl2m is the nuclear magneton and 
i^n = 3n/2, Kp = gp/2 — 1 with = —3.8263 and 
gp = 5.5856 being the intrinsic magnetic moments of pro¬ 
tons and neutrons. Interactions with the external mag¬ 
netic field are marked by the superscript (ext). This field 
is considered to be externally generated, and therefore 
there is no associated field equation and thus no other 
bosonic terms in the Lagrangian. 

The magnetic held breaks spherical symmetry for the 
Dirac and Klein-Gordon equations m- Only axial sym¬ 
metry is preserved. In addition, time-reversal symmetry 
is broken by the magnetic held leading to the appear¬ 
ance of time-odd mean helds and non-vanishing currents 
which induce space-like components of the vector mesons 
w and p, usually referred as nuclear magnetism [45H47j . 

The effects of the coupling of protons and neutrons to 
an external magnetic held can be classihed as follows: 
nucleon paramagnetism, caused by the interaction of the 
magnetic held with the neutron (proton) magnetic dipole 
moment. Since the gyromagnetic factor for neutrons 
(protons) is negative (positive), conhgurations with the 
spin anti-parallel (parallel) to the magnetic held are en¬ 
ergetically favored; and proton orbital magnetism, caused 
by the coupling of the orbital motion of protons with the 
magnetic held. It favors conhgurations where the pro¬ 
ton angular momentum projection is oriented along the 
direction of the external magnetic held. In general, it 
is thus expected that the magnetic held effects on the 
single-particle structure of nuclei are more pronounced 
for protons than for neutrons. 
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FIG. 1. (Color on-line) Binding energy BE predicted by DD- 
ME2 as a function of the strength of the magnetic held B 
for several typical nuclei appearing in the outer crust. A 
quadratic ht is also shown to indicate the average growing 
dependence of BE. 



From a more qualitative point of view, we show in Fig. 
I^the trends predicted by DD-ME2 in the binding energy 
BE, entering directly in the calculation of M{A,Z), as 
a function of the external magnetic held for several typ¬ 
ical nuclei thought to be present in the outer crust of 
a neutron star. From this hgure, one sees on average a 
parabolic increasing trend of the binding energy with the 
magnetic held -for guidance, htted parabolas are also 
showrQ We have checked that on average, the bind¬ 
ing energy of nuclei present in the outer crust does not 
increase by more than a 10 % when the more extreme 
magnetic helds {B ~ G) are taken into account. 

In the present work, we neglect the effect of pairing cor¬ 
relations on the nuclear binding energy, as it is expected 
that pairing effects will decrease with increasing magnetic 
held [44] . Indeed, the interaction of the nucleus with the 
external magnetic held removes all degeneracies in the 
single-particle spectrum, and tends to separate formerly 
degenerate levels with opposing signs of angular momen¬ 
tum projection (cf. Figs. 2 and 3 of Ref. [37]). Such 
single-particle energy splitting will produce a reduction 
of the neutron and proton pairing gaps with increasing 
magnetic helds and, eventually, their disappearance]^ 


^ For a given unit cell (i.e. fixed volume) and a uniform B, the 
classical magnetic field energy scales with 
^ In analogy, the nuclear cranking model has been very useful for 
describing some transitional nuclei that behave as a rotor, thus, 
displaying a collective rotation of frequency Urot- In this context, 
it can be shown that for large values of the rotation frequency of 
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C. Effects of different lattice configurations 


The lattice energy is not directly affected by the mag¬ 
netic held [48]. Nevertheless, some indirect effects on the 
lattice conhguration due to Coulomb screening of the ions 
embedded in the electron gas may arise. The calculation 
of the potential energy of the Coulomb lattic^ consists of 
divergent contributions that must cancel out as required 
by the overall charge neutrality of the system. Accurate 
calculations for the electron gas have been available for 
a long time and these results can be generalized to our 
case. Specihcally, in the absence of a magnetic held, it 
has been shown that the most energetically favorable con¬ 
hguration is a crystallization into a body-centered cubic 
lattice [50l|5T|. The lattice energy per baryon ei{A,Z) 
may be written as m, 

ei{A,Z,n) = -(1.81962)^^ = (13) 

where a is the lattice constant, e is the electron charge, 
Pp is in MeV and we have dehned the coefhcient 
Cbcc = 3.40665 X 10“^. The ei dependence on the 
density n or, equivalently, on the average baryonic 
Fermi momentum pp, enters via the lattice constant 
a. Similar calculations can be carried out evaluating 
diherent lattice conhgurations, like faced-centered cubic 
or simple cubic ones: ei dependencies on A, Z, n hold 
and only Cbcc has to be replaced by Cfcc or Csc [50] . 
The two latter coefficients are both smaller than Cbccj 
suggesting that a body-centered cubic lattice is the most 
favorable conhguration in the absence of a magnetic held. 


The body-centered cubic lattice pressure can be writ¬ 
ten as, 


Pi{A,Z) = -^ 

U ; gy 


Z,A 


2^/ 2 f^bcc • (14) 


In the absence of a strong magnetic held, the effect 
on the lattice conhguration due to the Coulomb screen¬ 
ing of an ion embedded in a uniform electron gas can be 
neglected as far as the composition of the outer crust is 
concerne In the presence of strongly quantiz¬ 

ing magnetic helds such as those found in magnetars, a 


about /icjcrit ~ 0.45 MeV where the pairing gap disappears m, 
the level degeneracy breaking also occurs. If one naively relates 
this rotation to be caused by an external magnetic field, one may 
find that cBcrit ~ o^crit^c/e = / \/^ ~ 5 X 10^ MeV^ 

which corresponds to about 10^^ G. 

^ We would like to note that the quantum zero-point motion of 
ions m has been neglected throughout this work. 

^ Indeed, in the absence of a magnetic field and considering the 
specific conditions present in the outer crust, both the elec¬ 
tron screening length Ae ~ ^ -within the mass-less 

Thomas-Fermi approximation El EH- and the ion-ion separa¬ 
tion rion-ion = ^ decrease with increasing 


recent work [^ has shown that the most favorable lattice 
configuration for nuclei may not be the body-centered cu¬ 
bic lattice. Intense magnetic fields cause an anisotropic 
screening of the Coulomb force by the electron gas lead¬ 
ing to Friedel oscillations in the ion-ion potential. Hence, 
for different values of the magnetic field, the authors of 
Ref. [S0] show that different Coulomb lattices such as 
faced-centered, hexagonal close-packed (hep) or body- 
centered cubic oriented along the magnetic field could 
emerge. We have estimated that, for our purposes, the 
energy differences do not essentially affect the predicted 
structure and composition of the outer crust. On the 
other hand, it has been shown that, under certain con¬ 
ditions m, interpenetrating cubic lattices formed by 
different ions is energetically favourable with respect to 
a bcc lattice of any other single ion (assuming a uni¬ 
form electron gas and in the absence of a magnetic field). 
Qualitatively, the same features are expected to remain 
when anisotropies of the background electron gas appear 
due to, for example, intense magnetic fields. Therefore, 
we have neglected the small energy correction due to such 
effects. 

Among the three terms in Eq. the lattice configu¬ 
ration -also considering screening, anisotropies or other 
secondary effects- do not significantly affect the crust 
composition. We verified that employing different lattice 
configurations or even neglecting the lattice for different 
magnetic values, the composition obtained is very similar 
with respect to the body-centered cubic case. In this re¬ 
gard the lattice configuration plays a very minor role and 
therefore from now we will simply treat nuclei as vertexes 
of a body-centered cubic lattice. 

For a fixed uniform magnetic field and pressure, the 
calculation of the Gibbs energy per particle can be now 
explicitly written as, |21] 


M{A,Z\B) 


+ (4 




MA,Z) 


g{A,Z-P,B) = 

^ \ ltf> 

(15) 

and one can search for the optimal nucleus, solving the 
set of equations that determine pe, ^max, rie, pFe{v)'. 


— F7.g ^I -t- 2^'aiax.R*^ 

PFciyY -I- TOg -I- 2vB^}j = pI 0 < F < 

^ 2^ 9uXe{l^) 


(16) 


Tie = 

{P = Pe + Pl{A, Z) 


electron densities. In particular, ~ 

for the whole outer crust. This considerations justify the no¬ 
screening approximation. 
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D. The neutron-drip transition point 


In this subsection we schematically study, using a 
toy model, how the pressure at the neutron-drip point 
changes in the presence of extreme magnetic fields {B > 
10^®G) which has both a direct effect on the pressure, 
and an indirect one through changes in the nuclear bind¬ 
ing energies -nuclei do not contribute to the pressure in 
our model but determine the electron chemical potential 
at the bottom layer of the outer crust. For more details 
on the effects on the magnetic field on the neutron-drip 
transition point, we address the reader to Refs. |22l I23j . 
In these works, where the magnetic field effect on nuclei 
was neglected, it was found that Pdrip increases linearly 
with extreme magnetic fields {B > 10^®G). 

To develope a simple, yet physical, model of the 
neutron-drip transition pressure Pdrip > we will neglect the 
small lattice contribution to the Gibbs free energy per 
baryon and, consistently also to the pressure. This as¬ 
sumption applied into the neutron-drip transition point 


allows us to write Eq.(15) as follows 


Ad' 


Nd BEd Zd 


(17) 


where the subscript’d’ denotes that the quantity should 
be evaluated at the neutron-drip transition. After some 
algebra, one finds, fj,e,d ~ nT-n — iTip + BEd/Zd where the 
neutron to proton mass difference can be neglected at 
this approximation level. 

Based on the recent results found in Ref. [23] and 
within our approximations, one may write the following 
expression for the dripping pressure. 


Pr 


drip 




Bi,m\ BE^ 






(18) 


that is valid only for extreme magnetic fields {B > 10^® 
G). Eq.(18) suggest that Pdrip changes linearly with the 
magnetic field and quadratically with the binding energy 
of the drip nucleus, which, at the same time, increases 
with increasing magnetic field (cf. Fig0. Specifically, 
we expect from this formula together with the average 
results shown in Figthat the increase of Pdrip between 
B = 10^® G and B = 10^^ G should be linear with P* 
and increase (roughly) one order of magnitude since the 
binding energy in nuclei is just barely affected on average 
by the strong magnetic field |23| . The situation we expect 
between B = 10^^ G and B = 10^® G is similar since 
the binding energy increases on average by a few % (no 
more than a ^ 10%, cf. Fig. 1^. Therefore, Pdrip should 
increase according to Eq. @1 ry one order of magnitude 
corrected by a small factor due to the increase in BE. 


III. RESULTS 

We present in this section the main results obtained 
from the numerical minimization of the Gibbs energy per 



FIG. 2. (Color on-line) Proton and neutron trends (a) and 
proton fraction y{P) = Z{P)/[N{P) + Z{P)] trends (b) pre¬ 
dicted by Duflo-Zuker model (blue solid line), by DD-ME2 
model (red dashed line) and by NL3 model (green dash-dotted 
lines), in absence of magnetic field. 


particle, g{A,Z\P,B), from the outermost part of the 
outer crust (ionization at a pressure PiorQ to the inner¬ 
most part (neutron-drip transition at a pressure PdrijQ. 
Beyond the neutron-drip transition, where the inner crust 
begins, the presence of the neutron gas requires the ad¬ 
dition of an extra term in the Gibbs energy 
We focus our study on the outer crust composition as a 
function of the pressure P(n), paying special attention 
to the effects produced by strong magnetic fields. To un¬ 
derstand the importance of the change in binding energy 
of nuclei induced by the magnetic field, results including 
(Pnuci = B) and excluding its effects (Pnuci = 0 G) on 
the composition will be presented alongside. The effect 
of the magnetic field on electrons is always taken into 
account. 

If nuclear binding energies are assumed not to be af¬ 
fected by the external magnetic field (Pnuci = OG), one 
expects to find the same nuclear composition than in the 
P = 0 G case with the only difference that the transitions 
from one nuclear species to another occur at higher pres¬ 
sures m US]- The latter effect is due to the electron 
interaction with the magnetic field: it reduces the Fermi 
energy /Ze and therefore delays, in terms of pressure, the 
appearance of a new nuclear species. It is important to 


® The ionization pressure Pjon corresponds, in good approxima- 
tion, to the free Fermi electron density at which the electron 
binding to the nucleus SF^eiec is not any more favorable. That 
is, BEeiec{Z) = At low pressure (density) the most fa¬ 

vorable nucleus has Z ~ 26 — 28, this implies a Pion ~ 10“^^ 
MeV fm-3. 

® The neutron-drip transition is determined by the condition 
g{Pdrip) = 'mn‘ it corresponds to the condition in which it is 
energetically favorable for the system to start dripping neutrons 
from the nucleus and form a neutron gas. 
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FIG. 3. (Color on-line) Compositions, i.e. Z{P) and N{P) 
trends, of the outer crust of a magnetar employing DD- 
ME2 model for three external magnetic fields B values: (a) 
B = 10^'‘G, (b) B = 10^® G, (c) B = 10^® G. Solid (blue) 
lines include effects of magnetic field on nuclei and electrons, 
i.e. Bnuci = B. Dashed (red) lines consider only magnetic 
effects on electron gas only, i.e. Bnuci = 0 G. Note that upper 
lines correspond to N{P) and lower lines to Z{P). 


FIG. 4. (Color on-line) Compositions, i.e. Z{P) and N{P) 
trends, of the outer crust of a magnetar employing NL3 model 
for three external magnetic fields B values: (a) B = lO^'^ G, 
(b) B — 10^® G, (c) B — 10^® G. Solid (blue) lines include ef¬ 
fects of magnetic field on nuclei and electrons, i.e. Bnuci = B.. 
Dashed (red) lines consider only magnetic effects on electron 
gas only, i.e. Bnuci = 0 G. Note that upper lines correspond 
to N(P) and lower lines to Z{P). 


note that for extremely large magnetic fields, of the or¬ 
der of i? ~ G, the effects of the electron gas can 

also change the composition [^. If one also accounts 
for the effects of a strong magnetic field on the nuclear 
binding energies (iJnuci = B ^ 0), further changes in the 
composition and discrepancies in the appearance of the 
subsequent nuclear species might be expected. As dis¬ 
cussed in Sec. |IIB[ as B increases, nuclei display larger 
binding energies on average (see Fig.[^. This induces the 
appearance of different nuclei populating the Coulomb 
lattice. 

To check the qualitative behavior of the NL3 and DD- 
ME2 models, in Fig. we compare their optimal pro¬ 
ton Z{P) and neutron N{P) numbers as a function of 


pressure in the absence of magnetic field with the pre¬ 
dictions of the Duflo-Zuker model [55H5T] . one of the 
most accurate mass models available in the literature 
(root mean square deviation of 400 keV). We notice that 
NL3 and DD-ME2 follow the main trends predicted by 
Duflo-Zukei[^ That is, two main changes in the neutron 
number to a two stable plateaus and, in coincidence, a 


^ DD-ME2 displays a dip in the neutron and proton number be¬ 
tween P ~ 10“^ — 10“® MeV fm“® and NL3 displays a gradual 
change in the neutron number as the pressure increases. Both 
are expected behaviors due to the fact that we neglect pairing 
correlations. As mentioned, we are interested on the effect of 































































clear change on the proton number with a subsequent 
decrease due to electron capture processes. Looking in 
more detail, the main differences arise in the transition 
of the neutron and proton numbers at pressures between 
10“® and 10“® MeV fm“^. Depending on the symme¬ 
try energy predicted by each model at an average nu¬ 
clear density, the transitions between the different neu¬ 
tron plateaus may appear at very different pressures [55] : 
the larger the symmetry energy, the earlier the transition 
takes place. Hence, DD-ME2 shows the larger symmetry 
energy at an average nuclear density, then Duflo-Zuker 
and finally NL3. Although the models used are accurate 
in the description of stable nuclei, the model dependence 
just seen in the results at B = 0 G will translate to 
the B 0 G cases and will allow us to estimate, to a 
sizable extent, the model dependence of our results. In 
the inset of the same figure. Fig. the proton fraction 
y{P) = Z{P)/[N{P) + Z{P)] is displayed, highlighting 
their good agreement. 

In Fig.j^and Fig.j^are displayed the compositions, i.e. 
Z{P) and N{P) trends, for three selected external mag¬ 
netic field values B, employing respectively the DD-MF2 
and NL3 models: (a) B = 10^"* G, (b) B = 10^® G, (c) 
B = 10^® G. Solid (blue) lines include effects of magnetic 
field on nuclei and electrons, i.e. Bnuci = B. Dashed 
(red) lines consider only magnetic effects on electron gas, 
i.e. Bnuci = 0 G — note that upper lines correspond to 
N{P) and lower lines to Z{P). 

We focus first on the DD-ME2 results for Bnuci = 0 
G -dashed (red) lines in Figj^, that is, only electrons 
are assumed to feel the magnetic field effects. We see in 
Figj^ that the changes in the neutron plateaus -from 
= 28 to = 50 and then to A^ = 82- are essentially 
unchanged with respect to the B = 0 G case. The same 
happens for protons, they follow the same trend. Re¬ 
garding the composition, it remains very similar to the 
B = 0 G case depicted in Fig. In the same panel, 
we then see, that turning on the effects of the magnetic 
field on nuclei -solid (blue) lines, the appearance of the 
low pressure neutron plateau is shifted towards higher 
pressures, leaving unchanged the high pressure one. In 
addition, the composition is slightly changed at the lower 
pressures shown in FigiB 

In Figj^, there is a clear shift of the neutron plateau 
A^ = 50 to higher pressures in both cases (Bnud = 0 
G and Bnud = B). The results on Bnud = 0 G are 
therefore in agreement with previous literature nansj. 
The predicted composition is similar in the two cases - 
except for a dip in the case in which Bnud is active- and 
similar to the B = 0 G case. Again, the main difference 
regarding the composition is on the low-pressure regime 
and due to the effect of the magnetic field on the binding 
energy of nuclei. 


very high magnetic fields on the outer crust where pairing effects 
are expected to be reduced. 

These changes are thought to be model dependent since we will 
see that NL3 predicts no change. 


For the strongest magnetic field, 10^® G, shown in 
Figj^, the situation is different. The structure between 
the Bnud = 0 G and Bnud = B cases is very similar and 
clearly different to the B = 0 G. One may say that just 
one nucleus is composing the outer crust for almost the 
whole range of pressures. Regarding the composition, it 
is now clear that including the magnetic field effects on 
nuclei will be very important for a precise understand¬ 
ing of the outer crust in such conditions. Specifically we 
found 4 §Zr 52 to be the nucleus in such a constant plateau 
if the magnetic field effects on nuclei are taken into ac¬ 
count. 

In order to estimate the model dependence of the re¬ 
sults we have just presented in Figj^ we show in Figj^ 
our results using a different model, the so called NL3. 
In Figj^, we show the results for B = 10^^ G, there 
are no essential differences between the Bnud = 0 G and 
^nud = B cases and with respect to the case in which 
the magnetic field effects on both electrons and nuclei are 
completely neglected (cf. Fig. for B = 0 G). 

Regarding the B = 10^® G predictions, shown in 
Figj^j, we see that the composition in the low pressure 
regime is different when compared to Figj^, for both 
.Bnud = 0 G and Bnud = B. The N = 50 plateau ap¬ 
pears shifted to higher pressures also in both cases. The 
case in which Bnud = 0 G -red (dashed) lines- follows 
the same trends as the Bnud = B G -blue (solid) lines- 
except for a region in which a dip in the composition is 
found. This region spans almost two orders of magnitude 
in pressure, from 3 x 10“^ MeV fm“® to about 10“® MeV 
fm“®. This behavior has been also seen in Figs. B and 
the details of which are model dependent. 

For the highest magnetic field considered here, shown 
in Figj^, we found the same type of behavior as in Figj^. 
That is, there is a constant plateau of neutron and pro¬ 
ton numbers for almost the whole range of pressures rel¬ 
evant for the outer crust. Gonsidering the effects of the 
external magnetic field on the binding energy of nuclei, 
turns out to be important. Specifically, we found 4 oZr 52 
again to be the nucleus in such a constant plateau for 
NL3 if the magnetic field effects on nuclei are taken into 
account. We have to emphasize that there is nothing spe¬ 
cial about this nucleus or in the fact that both models 
predict the same, being the change in single-particle level 
scheme the driving mechanism of the model dependence 
on the magnetic field, and considering that, for stable 
mid-mass nuclei, the level schemes are very similar for 
both DD-MF2 and NL3. In addition, it is expected that 
models with differing effective masses will yield a differ¬ 
ent nucleus sitting at the lattice. 

Summarizing, for field strengths of about G the 
inclusion of the effect of the magnetic field on nuclear 
binding energies may play a role (depending on the nu¬ 
clear mass model used, cf. Fig.[^ and Fig|^), and can¬ 
not be neglected for higher field strengths. If prominent 
magnetic fields, B = xo^®“^® G, are confirmed to exist in 
the the surface of neutron stars, one should expect some 
changes in the structure and composition of the outer 
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FIG. 5. (Color on-line) Neutron-drip transition pressures 
I^drip as a function of five external magnetic field B-. B ~ OG, 
B = 10^'‘G, B = lO^'^G, B = 10^® G, B = 10^'^ G and 
B = 10^® G. Circles (blue) refer to the results considering the 
effects of the magnetic field on both electrons and nuclei, i.e. 
Bnuci = B, while triangles (red) refers to calculations where 
the effects on nuclear binding energy have been neglected, i.e. 
I^nucl ~ 0. 


neutron-drip transition occurs, respectively, at a pressure 
one and two orders of magnitude higher than for weaker 
field strengths, where the Pdrip coincides with the results 
for zero magnetic field. The effect of the magnetic field on 
nuclei produces just a small change in Pdrip- Our results 
are in agreement with Refs. |221 [23l [25] since a change by 
an order of magnitude in B for the most extreme cases, 
induces a change of one order of magnitude in Pdrip, that 
is, (roughly) a linear change. It is not surprising, with 
Figj^and Eq.([I^ on hand, we could have expected the 
results shown in Figj^ changes on Pdrip induced by the 
effect of the magnetic field on nuclear binding energies 
are overshadowed by changes induced by the magnetic 
field on the electrons. 

In Fig. is displayed the equation of state P{n) as 
predicted by the DD-ME2 model, for six magnetic field 
values; for the NL3 model the behavior is similar. Panel 
(a) was obtained neglecting the effects of the magnetic 
field on nuclei, while on panel (b) they were fully taken 
into account. In the low density range, the magnetic 
field effect on the electron gas promotes a higher matter 
incompressibility, as signaled by the steep slope. As the 
pressure further increases ly —>■ and consequently 

the number of electronic levels which can be populated 
increases rapidly. Then, the magnetic field effects lose 
importance and the equation of state tends to approach 
the straight line associated to the B = 0 case. 


IV. CONCLUSIONS 


crust compared to the zero magnetic field case: 

• The effect of the magnetic field on electrons shifts 
to higher pressures the appearance of a new nuclear 
species, keeping almost unchanged the nuclei that 
populates the Coulomb lattice [HES]. 

• The magnetic field effect on nuclei favors heavier 
nuclei at lower pressures. As B increases, such 
heavier nuclei are also preferred up to higher pres¬ 
sures. In the most extreme case, and for the studied 
models, the whole outer crust is almost composed 
of 4 §Zr 52 , regardless of the nuclear model used to 
calculate the nuclear binding energies. 

• Extreme magnetic fields also favor in some cases the 
appearance of odd-mass nuclei in the outer crust, in 
contrast to what happens for lower magnetic fields 
where the pairing interaction has a more relevant 
role. 

In Fig. [^we display the evolution of the neutron-drip 
transition pressure Pdrip as a function of the external 
magnetic field. Circles (blue) refer to the results includ¬ 
ing the effects of the magnetic field on both electrons and 
nuclei, i.e. i?nuci = B, while triangles (red) indicate that 
the effects on nuclear binding energy have been neglected, 
i.e. Snuci = 0. For the strongest fields, the 


For the first time the composition of the outer crust of a 
cold non-accreting and strongly magnetized neutron star 
has been studied taking into account, in a self-consistent 
fashion, the effects of the magnetic field on the binding 
energies of the nuclear species present in the Coulomb 
lattice. Results both including and neglecting the effect 
of the magnetic field on nuclear binding energies have 
been presented in order to understand its impact. 

The required nuclear binding energies have been calcu¬ 
lated using the NL3 and DD-ME2 relativistic mean-field 
models with explicit and fully self-consistent couplings to 
the magnetic field. These two models have been proven 
to be able to reproduce a wide variety of experimental 
data, including binding energies, with a reasonable accu¬ 
racy. 

Important changes in the structure and composition 
of the outer crust in the limit of high-magnetic fields 
{B ~ G) have been found, even when neglecting 

its influence on nuclear binding energies. A shift to higher 
pressures, as the magnetic field increases, is observed in 
the transition from one nuclear species to another. The 
magnetic field effect on nuclear binding energies favors 
the appearance of heavier nuclei at low pressure, and, as 
the magnetic field increases, those heavier nuclei tend to 
be preferred up to greater pressures. In the most extreme 
case, and for the studied models, almost the whole outer 
crust is composed of 4 oZr 52 . Extreme magnetic fields also 
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FIG. 6. (Color on-line) Equations of state (n versus P) pre¬ 
dicted by the DD-ME2 model, for six different external mag¬ 
netic field values {B = OG, B = G, B = 10^® G, B = 
10^® G, B = 10^’’ G and B = 10^® G), including (Bnud = B) 
and neglecting (Bnuci = 0 G) the effects of magnetic fields on 
nuclei (respectively a and b). Curves are plotted up to the 
neutron-drip transition point. 


favor in some cases the appearance of odd-mass nuclei in 
the outer crust, in contrast to what happens for lower 
magnetic fields where the pairing interaction has a more 
relevant role. Furthermore, the neutron-drip transition 
pressure Pdrip for the highest magnetic fields considered 
{B > lO^^G) is increased from one and up to two orders 
of magnitude with respect to the no magnetic field case. 
The latter will impact the spatial extension of the outer 
crust. 

In summary, the inclusion of the effect of the magnetic 
field on nuclear binding energies may play a role in the 
determination of the outer crust composition, and thus 
in its mechanical and thermodynamical properties, for 
field strengths of about 10^^ G (depending on the nuclear 
mass model used), and cannot be neglected for high field 
strengths {B > 10^® G). 
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Appendix: numerical example 

A complete sequence of the outer crust composition, 
taken into account the effect of the magnetic field on nu¬ 
clei, is shown as an example in Table|T] Pmin (Pmax) is the 
minimum (maximum) pressure, in units of MeV fm“^, at 
which the given nucleus is present. The density Umax, 
expressed in units of fm“^, is the maximum baryonic 
density at which the given nucleus is present, according 
to Eq. (|^ and to charge neutrality. Pjon is the electronic 
ionization pressure, which represents the lower pressure 
limit of the outer crust and y is the proton fraction. The 
compositions correspond to the predictions of the DD- 
ME2 (left half side) and NL3 (right half side) models 
for a P = 10^® G case. It corresponds to the results in 
Fig|^ and Fig|^, respectively, depicted by using solid 
(blue) lines. 
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TABLE I. Composition of the outer crust of a magnetar with B — Bnuci = 10^® G. Pmin (Pmax) is the minimum (maximum) 
pressure, in units of MeVfm“®, at which the given nucleus is present. The density rimax, expressed in units of fm“®, is the 
maximum baryonic density at which the given nucleus is present, according to Eq. ^ and to charge neutrality hypotesis. Pion 
is the electronic ionization pressure, which represents the lower limit of the outer crust. 


NuclsUS -fmin -Pm ax 

^max 
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Nucleus Pmin Pmax 

^max 

y 

MeVfm-^ MeVfm-^ 
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fm-® 
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P. 
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0.464 

asSr 

1.44- 

10-" 

2.41 • 

10“® 

3.36- 

10"® 

0.448 

asSr 

1.72 • 

10-" 

2.00- 

10“® 

3.05- 

10"® 

0.432 

llRh 

2.42 • 

10"® 

2.45 • 

10“® 

3.44 • 

10"® 

0.431 

3?Rb 

2.01 • 

10“® 

2.18 • 

10“® 

3.22 • 

10"® 

0.425 

28Ni 

2.46 • 

10"® 

9.21 ■ 

10“® 

6.64- 

10"® 

0.412 

seKr 

2.19- 

10“® 

3.47- 

10“® 

4.09- 

10"® 

0.419 

alSe 

9.22 • 

10"® 

2.02 • 

10“® 

9.96- 

10"® 

0.409 

asBr 

3.48- 

10~® 

4.10 • 

10“® 

4.50- 

10"® 

0.412 

3oZn 

2.03- 

10"® 

2.12 • 

10“® 

1.06- 

10"® 

0.390 

34Se 

4.11 • 

10~® 

1.73- 

10“® 

9.23- 

10"® 

0.405 

'^®7n 

30^^ 

2.13 • 

10"® 

2.94- 

10“® 

2.23 • 

10"® 

0.385 


1.74- 

10“® 

2.36- 

10“® 

1.60- 

10"® 

0.390 

3oZn 

2.95- 

10"® 

4.09- 

10“® 

2.89- 

10"® 

0.380 

aiGa 

2.37- 

10“® 

2.73- 

10“® 

2.07- 

10"® 

0.383 

28Ni 

4.10 • 

10"® 

8.62- 

10“® 

5.65- 

10"® 

0.363 

3oZn 

2.74- 

10~® 

3.70- 

10“® 

2.73- 

10"® 

0.375 

2 lNi 

8.63- 

10"® 

1.44- 

lO”'* 

8.16- 

10"® 

0.354 

iNi 

3.71 • 

10~® 

1.44 • 

10“'‘ 

8.06- 

10"® 

0.359 

125rp 

43 1C 

1.45 • 

10-^ 

1.57- 

lO-'* 

9.61 • 

10"® 

0.344 

1227 

1.45 • 

10-4 

1.96- 

10-'‘ 

1.12 • 

10"^ 

0.328 

' 42 M 0 

1.58- 

10-^ 

1.74- 

lO”"* 

1.05- 

10-^ 

0.339 

121 v 
39 ^ 

1.97- 

10“^ 

2.26 • 

10“'‘ 

1.30- 

10"^ 

0.322 


1.75- 

10"^ 

1.90- 

lO-"* 

1.13- 

10"^ 

0.333 

^38Sr 

2.27- 

10“^ 

3.20- 

10-'‘ 

1.78- 

10"^ 

0.317 

1227 

1.91 • 

10"'^ 

2.70- 

lO”'* 

1.52 • 

10"'^ 

0.328 

“®Rb 

3.21 • 

10~^ 

3.35- 

10“'‘ 

1.86- 

10""^ 

0.311 

121 v 

39 ^ 

2.71 • 

10-^ 

3.08- 

10“^ 

1.70- 

10-^ 

0.322 

“sBr 

3.36- 

10~^ 

4.30- 

10-^ 

2.37- 

10""^ 

0.299 

^38Sr 

3.09- 

10-^ 

4.46- 

10“^ 

2.23 • 

10-^ 

0.317 

28Ni 

4.31 • 

IQ-'^ 

5.19- 

10-^ 

2.78- 

10""^ 

0.286 

“?Rb 

4.46. 

10"^ 

4.78- 

10“^ 

2.44 • 

lO"'^ 

0.311 



















12 


[13] N. V. Ardeljan, G. S. Bisnovatyi-Kogan, and S. G. Moi- 
seenko, Mon. Not. Roy. Astron. Soc. 359, 333 (2005). 

[14] D. Lai and S. L. Shapiro, Astrophys. J. 383, 745 (1991). 

[15] G. Y. Cardall, M. Prakash, and J. M. Lattimer, Astro¬ 
phys. J. 554, 322 (2001). 

[16] P. Haensel, A. Y. Potekhin, and D. G. Yakovlev, Neutron 
Stars 1: Equation of State and Structure (Springer, New 
York, 2007). 

[17] A.Y. Potekhin, and D.G. Yakovlev, Astron. Astrophys., 
314, 341 (1996), Astron. Astrophys., 327, 442 (1997). 

[18] A. Y. Potekhin, Astron. Astrophys. 351, 787 (1999). 

[19] J. Ventura and A. Y. Potekhin, Neutron Star Envelopes 
and Thermal Radiation from the Magnetic Surface, in 
Kouveliotou, G., Ventura, J., and van den Heuvel, 
E., eds., The Neutron Star - Black Hole Connection, 
Proceedings of the NATO Advanced Study Institute, 
Elounda, Grete, Greece, 7-18 June 1999, NATO Science 
Series G, vol. 567, pp. 393-414, (Kluwer Academic Pub¬ 
lishers, Dordrecht, Netherlands; Boston, U.S.A., 2001). 

[20] B. Carter, E. Chachoua, and N. Charnel, Gen. Rel. Grav. 
38, 83 (2006). 

[21] S. L. Shapiro and S. A. Teukolsky, Black Holes, White 
Dwarfs, and Neutron Stars: The Physics of Compact Ob¬ 
jects (Wiley, New York, 1983). 

[22] N. Charnel, A. F. Fantina, J. L. Zdunik, and P. Haensel, 
Phys. Rev. C 91, 055803 (2015). 

[23] N. Charnel et al, Phys. Rev. C 91, 065801 (2015). 

[24] J. M. Pearson, N. Charnel, S. Goriely and C. Ducoin, 
Phys. Rev. C 85, 065803 (2012). 

[25] N. Charnel et al, Phys. Rev. C 86, 055804 (2012). 

[26] G. Baym, C. Pethick, and P. Sutherland, Astrophys. J. 
170, 299 (1971). 

[27] S. B. Ruester, M. Hempel, and J. Schaffner-Bielich, Phys. 
Rev. C 73, 035804 (2006). 

[28] X. Roca-Maza and J. Piekarewicz, Phys. Rev. C 78, 
025807 (2008). 

[29] Pearson et al, Phys. Rev. C 83, 065810 (2011) 

[30] R. N. Wolf et al., Phys. Rev. Lett. 110, 041101 (2013) 

[31] S. Kreim et al, Int. J. Mass Spectrom., 63 349-350 (2013) 

[32] G. A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C 55, 
540 (1997). 

[33] J. Boguta, and A.R. Bodmer, Nucl. Phys. A 292, 413-428 
(1977). 

[34] W. Pannert, P. Ring, and J. Boguta, Phys. Rev. Lett. 
59, 2420 (1987). 


[35] G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, 
Phys. Rev. C 71, 024312 (2005). 

[36] K. Kiuchi, K. Kotake, Mon. Not. Roy. Astron. Soc. 385, 
1327 (2008). 

[37] D. Pena Arteaga, M. Grasso, E. Khan, and P. Ring, Phys. 
Rev. C 84, 045806 (2011). 

[38] 1. 1. Rabi, Z. Phys. 49, 507 (1928). 

[39] S. E. Agbemava, A. V. Afanasjev, D. Ray, and P. Ring, 
Phys. Rev. C 89, 054320 (2014). 

[40] Y. K. Gambhir, P. Ring, and A. Thimet, Annals Phys. 
198, 132 (1990). 

[41] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, and P. 
Ring, Phys. Rept. 409, 101 (2005). 

[42] J. D. Bjorken and S. D. Drell, Relativistic Quantum Me¬ 
chanics (McGraw-Hill, New York, 1964). 

[43] Y. R. Shimizu, J. D. Garrett, R. A. Broglia, M. Gallardo, 
E. Vigezzi, Rev. Mod. Phys 61 131 (1989). 

[44] C. Kittel, Quantum Theory of Solids (Wiley, New York, 
1963). 

[45] U. Hofmann and P. Ring, Phys. Lett. B 214, 307 (1988). 

[46] W. Koepf and P. Ring, Nucl. Phys. A 511, 279 (1990). 

[47] A. V. Afanasjev and H. Abusara, Phys. Rev. C 81, 
014309 (2010). 

[48] J. H. Van Vleck, The Theory of Electric and Magnetic 
Susceptibilities (Oxford University Press, London, 1932). 

[49] D. A. Baiko, Phys. Rev. E 80, 046405 (2009). 

[50] R. A. Coldwell-Horsfall and A. A. Maradudin, J. Math. 
Phys. 1, 395 (1960). 

[51] A. L. Fetter, J. D. Walecka, Quantum Theory of Many 
Particle Systems (McGraw-Hill, New York, 1971). 

[52] G. Watanabe and K. lida, Phys. Rev. C 68, 045801 
(2003). 

[53] C. J. Pethick, D. G. Ravenhall, Ann. Rev. Nucl. Part. 
Sci. 45, 429 (1995). 

[54] R. Sharma and S. Reddy, Phys. Rev. C 83, 025803 

( 2011 ). 

[55] D. A. Shalybkov and D. G. Yakovlev, Astrophysics 27, 
562 (1987). 

[56] P. F. Bedaque, S. Mahmoodifar, and S. Srimoyee, Phys. 
Rev. C 88, 055801 (2013). 

[57] C. J. Jog and R. A. Smith, Astrophys. J. 253, 839 (1982). 

[58] Nandi et al, Astrophys. J. 736, 156 (2011) 

[59] J. Duflo, Nucl. Phys. A 576, 29 (1994). 

[60] A. P. Zuker, Nucl. Phys. A5 576, 65 (1994). 

[61] J. Duflo and A. P. Zuker, Phys. Rev. C 52, 23 (1995). 



